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ABSTRACT 

Context. Giant radio halos (RH) are diffuse Mpc-scale synchrotron sources detected in a fraction of massive and merging galaxy 
clusters. An unbiased study of the statistical properties of RHs is crucial to constrain their origin and evolution. 

Aims. We aim at investigating the occurrence of RHs and its dependence on the cluster mass in a SZ-selected sample of galaxy 
clusters, which is as close as possible to be a mass-selected sample. Moreover, we analyse the connection between RHs and merging 
clusters. 

Methods. We select from the Planck SZ catalogue (Planck Collaboration XXIX 2014) clusters with M > 6x at z = 0.08-0.33 

and we search for the presence of RHs using the NVSS for z < 0.2 and the GMRT RH survey (GRHS, Venturi et al. 2007, 2008) and 
its extension (EGRHS, Kale et al. 2013, 2015) for 0.2 < z < 0.33. We use archival Chandra X-ray data to derive information on the 
clusters dynamical status. 

Results. We confirm that RH clusters are merging systems while the majority of clusters without RH are relaxed, thus supporting the 
idea that mergers play a fundamental role in the generation of RHs. We find evidence for an increase of the fraction of clusters with 
RHs with the cluster mass and this is in line with expectations derived on the basis of the turbulence re-acceleration scenario. Finally, 
we discuss the effect of the incompleteness of our sample on this result. 
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1. Introduction 

Clusters of galaxies are the largest and most massive bound sys¬ 
tems in the Universe. They form and grow at the intersection 
of cosmic filaments where matter and galaxies get together as 
a consequence of the gravitational collapse. Mergers between 
clusters of galaxies are among the most energetic events in the 
Universe as they release energies of ~ 10^^ - 10^^ erg in few 
Gyrs. Although most of this energy is dissipated to heat the in- 
tra cluster medium (ICM) up to temperature of ~ 10^ - 10^ °X, 
part of this energy is channelled into the acceleration of rela¬ 
tivistic particles and amplification of magnetic fields in the ICM 
(see e.g., Brunetti & Jones 2014, for a review). Diffuse Mpc- 
scale synchrotron radio emission observed in a growing number 
of galaxy clusters is the most direct and compelling evidence 
of this activity. Non-thermal radio emission from galaxy clus¬ 
ters is observed in the form of giant radio halos (RH), located 
at the cluster center with morphology similar to that of the X- 
ray emission, and radio relics, located at the cluster outskirts 
and characterised by elongated shapes (e.g., Feretti et al. 2012). 
The emerging theoretical picture is that radio relics trace shock 
waves propagating out of the cluster cores, whereas radio halos 
trace turbulent regions in clusters, where particles are trapped 
and re-accelerated during mergers (e.g., Brunetti & Jones 2014). 
The comparison between thermal and non-thermal properties 
of galaxy clusters provides important information on the com¬ 
plex mechanisms that generate the observed radio emission. 


According to models based on turbulent acceleration, the forma¬ 
tion history of RHs depends on the cluster merging rate through¬ 
out cosmic epochs and on the mass of the hosting clusters, which 
ultimately sets the energy budget available for the acceleration of 
relativistic particles. In their simplest form, these models predict 
a steepening in the spectra of RHs at a frequency which di¬ 
rectly depends on the energetics of the merger (i.e. on the cluster 
mass). Therefore a key expectation is that typical RHs should 
preferentially be found in massive objects undergoing energetic 
merging events, whereas they should be rarer in less massive 
merging-systems and absent in relaxed clusters (e.g., Cassano & 
Brunetti 2005). Smaller systems undergoing less energetic merg¬ 
ing events are expected to produce RHs with increasingly steep 
spectra (lower v^) which become underluminous at higher fre¬ 
quencies. This implies the existence of RHs with ultra-steep ra¬ 
dio spectra (USSRH, a > 1.5, with /(v) oc y~^) that should be¬ 
come better visible at low radio frequency (Cassano et al. 2006; 
Brunetti et al. 2008, Dallacasa et al. 2009). 

A first statistical measurement of the occurrence of giant 
RHs in galaxy clusters has been obtained through the “GMRT 
RH Survey” (Venturi et al. 2007, 2008; GRHS hereafter) and 
its extension (the EGRHS, Kale et al. 2013, 2015). This survey 
is restricted to clusters in the redshift range 0.2 - 0.4. It con¬ 
firmed that RHs are hosted in only ~ 20 - 30% of X-ray lumi¬ 
nous (Lx(0.1 - 2.4 keV) > 5 x 10"^"^ erg/s) clusters and found 
that clusters branch into two populations: RHs trace a corre- 
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lation between P\a and Lx, whereas radio-undetected clusters 
(upper limits) lie about 1 order of magnitude below the corre¬ 
lation {e.g., Brunetti et al. 2007). Importantly, this bimodal split 
can be traced to clusters dynamics: RHs are always associated 
to merging system^ while clusters without RH are typically re¬ 
laxed {e.g., Cassano et al. 2010). 

The recent advent of cluster surveys via the SZ effect 
(i.e., with the Planck satellite) has enabled the construction of 
unbiased cluster samples that are almost mass-selected, due to 
the tight relation between the total SZ signal, Fsoo, when in¬ 
tegrated within R 5 O (0 and the cluster mass, M 500 (Motl et al. 
2005; Nagai 2006). The fraction of clusters with RHs appears 
larger in SZ-selected cluster samples with respect to that derived 
from X-ray samples (Sommer & Basu 2014). Earlier studies 
were unable to observe a bimodal behaviour of clusters with RH 
and radio-undetected systems in the radio-SZ properties (Basu 
2012). However, thanks to the improved statistics, more recently 
Cassano et al. (2013) demonstrated the presence of a bimodal 
split also in the radio-SZ diagram, for F 500 > 6 x 10“^ Mpc^, and 
confirmed that such split is tightly connected with the dynamical 
properties of the hosting clusters. This result provides strong ev¬ 
idence, complementary to X-ray studies, that mergers play a key 
role in the formation of RHs. However, the relatively low mass 
completeness (~ 50%) of the Cassano et al. (2013) sample did 
not allow to measure the occurrence of RHs and in particular to 
study such occurrence as a function of the cluster mass. 

In order to provide an unbiased measure of the fraction 
of clusters hosting RHs and of its dependence on the clus¬ 
ter mass, we selected from the PSZ catalogue clusters with 
M > 6 X 10^“^ Mq in the redshift range z ^ 0.08 - 0.33. In this 
way we obtained a sample of 75 clusters with mass complete¬ 
ness > 80%, 57 of them have available radio information. Here 
we report on the statistical analysis of these 57 clusters that con¬ 
stitute a sample with mass completeness ~ 63% (see Sect.|^and 
|7]). The addition of the remaining 18 clusters without radio in¬ 
formation will allow to achieve a completeness in mass > 80%. 
Deep JVLA and GMRT observations of these 18 clusters are in 
progress and results will be presented in a follow-up paper. 

We also used the Chandra X-ray data, available for most of the 
clusters in the sample, to investigate their dynamical status and 
the connection with the radio properties. 

In Sect. 2 we describe the selection of the cluster sample; in 
Sect. 3 we report on the analysis of NVSS data of low-z clusters; 
in Sect. 4 we derive the cluster dynamical status. In Sect. 5 we 
derive the occurrence of clusters with giant radio halos and in 
Sect. 6 we investigate the RH-merger connection. In Sect. 7 we 
discuss the effect of the sample completeness on the results and 
we report our conclusions in Sect. 8 . 

2. Cluster sample selection 

We used the Planck SZ cluster catalogue (PSZ, Planck 
Collaboration XXIX 2014) to select a sample of massive galaxy 
clusters. This catalogue consists of 1227 objects derived from SZ 
effect detections using the first 15.5 months of Planck satellite 
observations. It contains 861 confirmed clusters and 366 cluster 
candidates. To date the Planck sample is the largest SZ-selected 
cluster sample (six times the size of the Planck Early SZ, Planck 
Collaboration 2011) and the deepest all-sky catalogue. It spans 


^ A possible outlier is the RH recently discovered in the cool-core 
cluster CL1821+643 (Bonafede et al. 2014). 

^ R 500 is the radius corresponding to a total density contrast 500pc(z), 
Pc{z) is the critical density. 



Fig. 1. Average mass limit computed from the average noise over the 
sky for the PSZ catalogue. The dotted, dashed and solid lines show the 
Planck mass limit at 80, 50 and 20 % completeness respectively. The 
rectangles show the regions from where we extracted our sample: the 
red rectangle for the low redshift sample, the cyan retangle for the high 
redshift one. Adapted from Planck Collaboration XXIX 2014. 


the broadest cluster mass range from 0.1 to 1.6 x 10 ^^ Mq, with 
redshift up to about one. 

Prom the PSZ catalogue we selected clusters with M 500 ^ 
6 X 10^'^M(£]and redshift 0.08 < z < 0.33. To maximize the 
radio coverage we adopted a declination limit 6 > -31°and 
\b\ > 20 °(|Z 7 | is the galactic latitude) for clusters at z > 0 . 2 , that 
coincides with that of the CMRT radio surveys. On the other 
hand at lower redshift we adopted 6 > -40° to ensure a follow 
up from the NVSS radio survey (Condon et al. 1998). 

Among the 54 clusters at z > 0.2, 34 belong to the EGRHS 
and thus have deep radio observations, additional 2 clusters have 
literature information (namely PSZl G205.07-62.94 and PSZl 
G171.96-40.64). 

Por targets in the redshift range 0.08-0.2 we collected data from 
the literature (14 clusters) and analysed data from the NVSS ra¬ 
dio survey (Condon et al. 1998) for the remaining 7 clusters (see 
Sect.[^. 

In Pig[^we show the M 500 -z distribution of the Planck clus¬ 
ters detected over 83.7% of the sky, together with the Planck 
mass limit corresponding to the 80, 50 and 20% completeness of 
the catalogue (Planck Collaboration XXIX 2014). The red and 
blue boxes show the regions from where we selected our cluster 
sample: the low-z sample and the high-z sample have a mass¬ 
completeness of -90% and 80%, respectively. The sample with 
radio information consists of 57 clusters (21 at z = 0.08 - 0.2 
and 36 at z = 0.2 - 0.33) with a completeness in mass of ~ 90% 
at low redshift and - 53% (0.8 x = 0.53) at higher redshift. 
The 57 clusters and their properties are listed in Tab.[2 

3. The low-z sample and the NVSS data analysis 

We use the NRAO VLA Sky Survey (NVSS, Condon et al. 1998) 
to investigate the presence of cluster-scale diffuse emission in 
the 7 clusters of the low-z sample missing radio information in 
the literature. The NVSS is a radio survey performed at 1.4 GHz 
with the Very Large Array (VLA) in D and DnC configuration. 


^ The values of M 500 in the PSZ catalogue are obtained from F 500 as 
described in Sect.7.2.2 in Planck Collaboration XXIX 2014. 
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Table 1. Total sample clusters properties 


cluster name 

RA 

Dec 

z 

Msoo 

(10“*Me) 

Radio info 

X-ray info 

A1437 

12 00 22.3 

+03 20 33.9 

0.134 

5.69 

noRH* 

MV 

A2345 

21 27 06.8 

-12 07 56.0 

0.176 

5.71 

Relics'* 

MV 

A2104 

15 40 08.2 

-03 18 23.0 

0.153 

5.91 

noRH* 

MV 

Zwcl 2120.1-^2256 

21 22 27.1 

+23 11 50.3 

0.143 

5.91 

no RH‘ 

MV 

RXC J0616.3-2156 

06 16 22.8 

-21 56 43.4 

0.171 

5.93 

noRH* 

MV 

A1413 

11 55 18.9 

+23 24 31.0 

0.143 

5.98 

MH^ 

RV 

A1576 

12 37 59.0 

+63 11 26.0 

0.302 

5.98 

UL^ 

R” 

A2697 

00 03 11.8 

-06 05 10.0 

0.232 

6.01 

UL^ 

R^ 

Z5247 

12 33 56.1 

+09 50 28.0 

0.229 

6.04 

RH’ 

MV 

Zwcl 0104.9-^5350 

01 07 54.0 

+54 06 00.0 

0.107 

6.06 

RH« 

- 

RXC 10142.0-^2131 

01 42 02.6 

+21 31 19.0 

0.280 

6.07 

UL^ 

R” 

A1423 

11 57 22.5 

-h33 39 18.0 

0.214 

6.09 

UL’ 

R” 

RXC J1314.4-2515 

13 14 28.0 

-25 15 41.0 

0.244 

6.15 

RH‘ 

M^ 

A2537 

23 08 23.2 

-02 11 31.0 

0.297 

6.17 

UL’ 

R” 

A68 

00 37 05.3 

+09 09 11.0 

0.255 

6.19 

UL’ 

MV 

A1682 

13 06 49.7 

+46 32 59.0 

0.226 

6.20 

RH’ 

M” 

A1132 

10 58 19.6 

+56 46 56.0 

0.134 

6.23 

noRH’ 

M V 

RXJ1720.1-^2638 

17 20 10.1 

+26 31 29.5 

0.164 

6.34 

MH** 

RV 

A781 

09 20 23.2 

+30 26 15.0 

0.295 

6.36 

UL’ 

M” 

A2218 

16 35 51.6 

+66 12 39.0 

0.171 

6.41 

RH’ 

> 

oo 

A3411 

08 41 55.6 

-17 29 35.7 

0.169 

6.48 

RH‘® 

M V 

Zwcl 0634.1-^4750 

06 38 02.5 

+41 41 23.8 

0.174 

6.52 

suspect* 

MV 

A3888 

22 34 26.8 

-37 44 19.1 

0.151 

6.67 

suspect* 

M?^.29 

A3088 

03 07 04.1 

-28 40 14.0 

0.254 

6.71 

UL^ 

R” 

A2667 

23 51 40.7 

-26 05 01.0 

0.226 

6.81 

UL^ 

R” 

A521 

04 54 09.1 

-10 14 19.0 

0.248 

6.91 

RHii.t/s 

M” 

A2631 

23 37 40.6 

+00 16 36.0 

0.278 

6.97 

UL’ 

M” 

A1914 

14 26 03.0 

+31 49 32.0 

0.171 

6.97 

RH*’ 

MV 

RXC J1504.1-0248 

15 04 07.7 

-02 48 18.0 

0.215 

6.98 

MH*’ 

R” 

A520 

04 54 19.0 

+02 56 49.0 

0.203 

7.06 

RH*-* 

M” 

A478 

04 13 20.7 

+ 10 2S 35.0 

0.088 

7.06 

MH*’ 

RV 

A773 

09 17 59.4 

+51 42 23.0 

0.217 

7.08 

RH*-* 

M” 

A1351 

11 42 30.8 

-h58 32 20.0 

0.322 

7.14 

RH*’ 

M V 

A115 

00 55 59.5 

+26 19 14.0 

0.197 

7.21 

Relic*'* 

MV 

A1451 

12 03 16.2 

-21 32 12.7 

0.199 

7.32 

suspect* 

M^ 

PSZl G205.07-62.94 

02 46 27.5 

-20 32 5.29 

0.310 

7.37 

no RH^ 

M^ 

A2261 

17 22 17.1 

+32 08 02.0 

0.224 

7.39 

UL’ 

R” 

RXCJ2003.5-2323 

20 03 30.4 

-23 23 05.0 

0.317 

7.48 

RH* 

M” 

A2552 

23 11 26.9 

+03 35 19.0 

0.300 

7.53 

RH?’ 

R? V 

A3444 

10 23 50.8 

-27 15 31.0 

0.254 

7.62 

MH’ 

RV 

S780 

14 59 29.3 

-18 11 13.0 

0.236 

7.71 

MH’ 

R” 

A2204 

16 32 45.7 

+05 34 43.0 

0.151 

7.96 

MH*’ 

RV 

A1758a 

13 32 32.1 

+50 30 37.0 

0.280 

7.99 

RH*’ 

M” 

A209 

01 31 53.0 

-13 36 34.0 

0.206 

8.17 

RH* 

M” 

A665 

08 30 45.2 

+65 52 55.0 

0.182 

8.23 

RH’ 

MV 

A1763 

13 35 17.2 

+40 59 58.0 

0.228 

8.29 

noRH’ 

MV 

RXC J1514.9-1523 

15 14 58.0 

-15 23 10.0 

0.223 

8.34 

RH*’" 

MV 

A1835 

14 01 02.3 

+02 52 48.0 

0.253 

8.46 

MH*** 

RV 

A2142 

15 58 16.1 

+21 13 29.0 

0.089 

8.81 

RH’-* 

M” V 

A1689 

13 11 29.5 

-01 20 17.0 

0.183 

8.86 

RH’o 

M” V 

A1300 

11 31 56.3 

-19 55 37.0 

0.308 

8.83 

RH’*" 

M” 

A2390 

21 53 34.6 

+ \1 40 11.0 

0.234 

9.48 

MH*’ 

R” 

A2744 

00 14 18.8 

-30 23 00.0 

0.307 

9.56 

RH*-* 

M” 

A2219 

16 40 21.1 

-^46 41 16.0 

0.228 

11.01 

RH*’ 

M” 

PSZl G171.96-40.64 

03 12 57.4 

+0S 22 10 

0.270 

11.13 

RH”’" 

M^ 

A697 

08 42 53.3 

+36 20 12.0 

0.282 

11.48 

RH”’*” 

M”.’* 

A2163 

16 15 46.9 

-06 08 45.0 

0.203 

16.44 

RH” 

M” 


RH = Radio Halo, MH = Mini-Halo, UL = Upper Limit, M=merger, R= relaxed. ^ Venturi et al. (2007), ^ Venturi et al. (2008), ^ Giovannini & 
Feretti (2000), Bonafede et al. (2009), ^ Govoni et al. (2009), ^ Kale et al. (2013), ^ Kale et al. (2015), ^ van Weeren et al. (2011), ^Giacintucci 
et al. (2014b), ^^van Weeren et al. (2013), ^^Brunetti et al. (2008), ^^Bacchi et al. (2003), ^^Giacintucci et al. (2011a), ^"^Govoni et al. (2001), 
Giacintucci et al. (2014a), Giacintucci et al. (2009), ^^Giovannini et al. (2006), ^^Giacintucci et al. (2011b), ^^Murgia et al. (2009), ^^Vacca et 
al. (2011), Venturi et al., 2013, ^^Giacintucci et al. (2013), ^^Feretti et al. (2001), ^"^Farnsworth et al. (2013), ^^Cassano et al. (2010), ^^Cassano 
et al. (2013), ^^Andersson & Madejski (2004), ^^Pratt et al. (2005), ^^WeiBman et al. (2013a), ^^Owers et al. (2011), ^^Girardi et al. (2006), 
^^Macario et al., 2010, Ultra Steep Spectrum RH, ^ candidate USSRH, ^Ferrari et al. (private communication), ^ visual inspection of 
XMM-Newton image, *NVSS data analysed in this paper, ^J X-ray Chandra data analysed in this paper. 
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It covers the sky north of ^ = -40°, it has an angular resolution 
of 45” and a surface brightness rms of ~ 0.45 mJy/beam. 

The low-z sample includes clusters with 0.08 < z < 0.2 
and M 500 > 5.7 X 10^"^ Mq. Radio interferometers suffer from 
the lack of sampling at short baselines, resulting in decreased 
sensitivity to emission on large spatial scales. For this reason we 
adopted a lower redshift limit of z > 0.08. Indeed, Farnsworth et 
al. (2013) showed that on scales >11 arcmin, that correspond to 
a 1 Mpc halo at z ~ 0.08, less than 50% of the total flux density 
is recovered with a NVSS snapshot observation. 

The upper redshift limit (z < 0.2) and the minimum mass are 
set by the angular resolution and sensitivity of the NVSS. The 
NVSS beam of 45" corresponds to ~ 150 kpc at z = 0.2 and this 
does not allow to separate discrete sources from residual diffuse 
emission at higher redshift. 

From Eq. 9 in Cassano et al. (2012), adopting the sensitivity 
and resolution of the NVSS, we derived the minimum F 1.4 of 
a detectable RH. The minimum mass M 500 = 5.7 x has 

been then derived assuming the P\a- ^500 correlation (Eq. 14 
in Cassano et al. 2013). 

With these selection criteria, the low-z sample is made of 
21 clusters. For 14 of these clusters we found information in the 
literature on the presence/absence of cluster-scale radio emission 
that are based on pointed VLA/WSRT observations. 

3.1. NVSS data analysis 

Here we describe the NVSS data analysis carried out to inves¬ 
tigate the presence of diffuse radio emission in the 7 clusters 
(marked with * in Tab. [T]) which lack literature radio informa¬ 
tion. To improve the quality of the radio images, i.e., to lower 
the rms noise and reduce the contribution of noise pattern, we 
reprocessed the NVSS fields of these 7 clusters. Data were anal¬ 
ysed using the NRAO Astronomical Image Processing System 
(AIPS). We calibrated the NVSS dataset and we obtained the im¬ 
ages of the pointings containing the cluster, then we combined 
them with the task FLATN. This procedure, known as the mosaic 
technique, is fundamental especially when the cluster falls at the 
border of the primary beam, because the signal to noise ratio de¬ 
creases with the distance from the pointing position. For the 7 re¬ 
processed clusters we reached an average rms^^ 0.25 mJy/beam, 
which is ~ 2 times better than the nominal NVSS noise. 

None of these clusters show clear diffuse cluster-scale ra¬ 
dio emission, however, we further investigated the possible pres¬ 
ence of residual emission in the central regions of these clusters. 
Specifically, we selected on each map a 1 Mpc sized circle cen¬ 
tred on the centroid of the cluster X-ray emission. With the task 
BLANK, we masked the discrete sources in the cluster that show 
contours at least at the 6cr level, then we measured the residual 
diffuse flux density (RDF, hereafter) in the circular region. We 
compared the RDF with the flux densities measured in other ar¬ 
eas of the same size taken around the cluster (3 for each cluster), 
i.e., “control fields”. In order to make a consistent comparison 
we normalized both the RDF and the control field flux densities 
(CFF, hereafter) to the number of pixels enclosed in a circle of 
1 Mpc diameter after masking the discrete sources. An exam¬ 
ple of this procedure is reported in Fig. applied to the case 
of Zwcl0634.1-1-4750. We stress that few NVSS beams corre¬ 
spond to 100-300 kpc. Consequently the use of the task BLANK 
in the case of relatively bright central sources is expected to re¬ 
move also diffuse emission on these scales. This is particularly 
problematic for the case of mini-halos that however are not the 
central focus of this paper. 


To test the reliability of this procedure, we also applied it to 
3 known RH clusters (A3411, A2218, Zwcl 0104.9-r5350) and 
to the mini-halo in RXJ1720.1-r2638 that belong to our sample, 
but have information from the literature. 

In Fig. ^ (left panel) we compare the average value of the 
CFF with the RDF for each cluster, whereas the offset between 
the RDF and the average CFF are reported in Fig. (right 
panel). Four clusters (A3888, Zwcl0634.1-h4750, A1451 and 
A2104) and the clusters with already known RHs (A2218, Zwcl 
0104.9-r5350 and A3411), show an excess at > 2cr level. We 
consider this as the threshold level to identify clusters with the 
possible presence of a RH. We note that all the clusters show 
a positive offset (at least -few mJy) between the RDF and the 
CFF (Fig. right panel). This however is likely due to residual 
contamination from faint cluster radio galaxies that are below 
the NVSS detection limit, rather than to diffuse flux on cluster 
scales (e.g., Farnsworth et al. 2013). For the sake of complete¬ 
ness, in Fig. (right panel, black asterisks) we also show the 
expected level of RH emission according to the Pi a - M 500 cor¬ 
relation (Cassano et al. 2013, Eq. 4.10). We note that the mini¬ 
halo in RXJ1720.1-r2638 does not result as an excess of diffuse 
emission. As explained above, this is due to the fact that a large 
fraction of the diffuse emission associated to the mini-halo is 
masked with the central bright radio galaxy. 

Although with our procedure based on NVSS we can iden¬ 
tify cases with suspect diffuse radio emission, we cannot confirm 
the presence of RHs in A3888, Zwcl0634.1-1-4750, A1451 and 
A2104. Deeper observations at low resolution (e.g., VLA array 
C or D), in order to have a good sensitivity to the diffuse Mpc 
scale emission, and at high resolution (e.g., VLA array A or B), 
to make an accurate subtraction of the individual sources from 
the u-v data, are necessary. 

4. Cluster dynamical status 

In this Section we report on the analysis of the dynamical prop¬ 
erties of clusters using Chandra X-ray data. A high fraction of 
clusters of the sample (50 out of 57) has Chandra archival data. 
24 of them already have dynamical information in the literature 
(Cassano et al. 2010, Cassano et al. 2013, see Tab. [T]). We 
produced the X-ray images of the remaining 26 clusters (marked 
with V in Tab. in the 0.5-2 keV band using CIAO 4.5 (with 
calibration files from CALDB 4.5.8). We adopted an automatic 
algorithm for the identification of point sources which were 
then removed from images. Each image was then normalized 
for the exposure map of the observation, which provides the 
effective exposure time as a function of the sky position exposed 
on the CCD. In our analysis we did not correct for the back¬ 
ground emission to treat the exposure-corrected images without 
introducing negative values in correspondence of pixels with 
zero counts. This procedure is sufficiently safe both because we 
are dealing with integrated quantities and since inside ^^^=500 
kpc (see below) the images are largely dominated by the signal 
associated to the cluster emission. Typically, using the back¬ 
ground estimates provided in Tab. 6.7 at http://asc.harvard.edu 
<http: //asc .harvard. edu/>/proposer/POG/html/chap 6 .html, we 
found that -95% (both in ACIS I and ACIS S) of the total 
counts in the 0.5-2 Kev band are from the cluster. 

Following Cassano et al. (2010, 2013), we studied the cluster 
substructures on the RH scale analysing the surface brightness 
inside an aperture radius Rap = 500 kpc, since we are inter¬ 
ested in the cluster dynamical properties on the scales where the 
energy is most likely dissipated. We used three main methods: 
the power ratios (e.g., Buote & Tsai 1995; Jeltema et al. 2005; 
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Fig. 2. Zwcl0634.1+4750 NVSS map, the contour levels are 0.66 x (-1,1,2,4,6,8,10,16,32,64) mJy The Icr level is 0.22 mJy In both 
panels the region where we extracted the flux densities are shown (solid circles). The central region has a diameter of 1 Mpc and is centred on 
the centroid of the X-ray emission; the other 3 areas are the so-called control flelds. With the task BLANK we masked the discrete sources in the 
central region and the one that falls in the lower control held (left panel). 



0.06 


:o.o4 


0.02 


control fields average flux [Jy] 



Fig. 3. Left panel. Diffuse flux density vs. control flelds average flux density. Symbols: colored dots are clusters with reprocessed NVSS datasets; 
black filled dots are known RH clusters and the black open dot is the mini-halo. The red line is the 1:1 line. Right panel Excess of diffuse flux with 
respect to the average control flelds flux densities. Black asterisks represent the expected radio power of the RH on the basis of the P 1.4 - M 500 
correlation (Cassano et al. 2013). 


Ventimiglia et al. 2008; Bohringer et al. 2010), the emission cen¬ 
troid shift {e.g., Mohr et al. 1993; Poole et al. 2006; O’Hara et al. 
2006; Ventimiglia et al. 2008; Maughan et al. 2008; Bohringer 
et al. 2010), and the surface brightness concentration parameter 
{e.g., Santos et al. 2008). 

The power ratio represents the multipole decomposition of 
the two-dimensional mass distribution inside a circular aperture 
Rap, centred on the cluster X-ray centroid. The power ratio can 
be defined as: 


Pq — [^0 


( 1 ) 


where gq is the total intensity inside the aperture radius: gq = 
S (< Rap), S (x) is the X-ray surface brightness, and 


Pm 


1 


(^m + ^m) 


( 2 ) 


where the moments and bm are given hy: 

am(R) =1 S {x'){R')cos{m<p')(f x' (3) 

■jR'<Rap 

and 

b„{R)= f S{x'){R’)sin{m<p’)d^x’ (4) 

\yR'<Rap 


we will only make use of the P 3 IP 0 parameter that is related to 
the presence of multiple peaks in the X-ray distribution provid¬ 
ing a clear substructure measure (Buote 2001, Bohringer et al. 
2010 ). 

The centroid shift, w is defined as the standard deviation of 
the projected separation between the peak and the centroid in 
unit of Rap and it is computed in a series of circular apertures 
centered on the cluster X-ray peak (e.g., Poole et al. 2006): 


w = 


N 


^ - {A}f 


1/2 


1 




(5) 
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Ai is the distance between the X-ray peak and the centroid of the 
ith aperture. 

Following Santos et al. 2008 we define the concentration pa¬ 
rameter as the ratio between the peak and the ambient surface 
brightness: 

_ S(r< 100 kpc) 

S(< 500 kpc) ^ 

The concentration parameter allows to distinguish clusters with 
compact core (not disrupted by recent mergers) from clusters 
with a spread distribution of the gas in the core. 

Basically, high values of P 3 IP 0 and w indicate a dynamically 
disturbed system, while high values of c stand for highly relaxed 
systems. 


5. Occurrence of Radio Halos 

The aim of this Section is to derive the occurrence of RHs as a 
function of the mass of the hosting clusters. Among the sample 
of 57 clusters with radio information, 24 host RHs and 4 show 
residual emission in a Mpc-scale regi on t hat is a possible indi¬ 
cation for the presence of a RH (Sect. |33- We split this sample 
into two mass bins and derived the fraction of clusters with RH, 
/rh, in the low mass bin (LM, M < Mum) and in the high mass 
bin (HM, M > Mum) for different values of the limiting mass, 
Miim (as detailed below). In general we found that /rh is lower in 
the LM bins (/rh ~ 20-30%) while it is higher (/rh ~ 60-80%) 
in the HM bins (Fig.[^. 

This difference is systematic and thus we attempted to iden¬ 
tify the value of that provides the most significant jump be¬ 
tween low and high-mass clusters. We performed Monte Carlo 
simulations considering both the cases in which the four objects 
in the low-z sample with suspect diffuse emission are included 
(/) as non RH clusters and (ii) as RH clusters. Considering the 
case (i), we randomly assigned 24 RHs among the 57 clusters 
of the sample and obtained the distributions of RHs in the two 
mass bins (after 10^ trials), expected in the case that RHs were 
distributed independently of the cluster mass. We consider 5 dif¬ 
ferent values of the transition mass between the two bins, specif¬ 
ically Miim = ( 6 ,7, 8 ,9,10) X 10^^ Mq. An example of the ex¬ 
pected distribution of the number of RH in the HM bin is shown 
in Fig.|^for the case Mum = 8 x 10^^ Mq. Each distribution can 
be well fitted by a gaussian function. The results of the Monte 
Carlo simulations are reported in Tab. for both the cases (i) 
(upper panel) and (ii) (lower panel). Specifically we report the 
number of clusters (Ndusters)^ the number of RHs (Nrh) and the 
fraction of clusters hosting RHs (fRu) in the two mass bins for 
each value of Mu^. In Tab. we also report the significance of 
our result in unit of cr, Z = (Nrh - p)Icr (where p is the gaussian 
median value), and the most likely value of /rh, fn=p IN duster, in 
the HM bin (very similar results are obtained for the LM bin). 

Fig. shows the observed fraction of RHs (dots) together 
with the results of the Monte Carlo simulations (shadowed re¬ 
gions) in the HM bin (left panel) and in the LM bin (right 
panel). We report the measured fraction of cluster with RHs 
and the results of the Monte Carlo analysis in both cases (i) 
(red and black dots and shadowed regions) and (ii) (green dots 
and shadowed regions). For a more clear visualization in Fig.[^ 
for the case (ii) we only show the results obtained by assuming 
Miim = (7, 8 ,9) X 10^"^Mq. Fig. shows that in the HM bin the 
observed fRu is always grater than that predicted by the Monte 
Carlo simulations, on the contrary in the LM bin the observed 



Npjj HM bin 

Fig. 4. Distribution of the number of RHs in the HM bin (M > Mum - 
8 X 10^"^Mq) after 10^ Monte Carlo trials. The red pont represents the 
observed number of RHs in the HM bin. 

fRH is always lower than that predicted by the Monte Carlo anal¬ 
ysis. This suggests the existence of a systematic drop of fRu in 
low mass systems. 

In both cases ( i) and ( ii), we found that the value of Mum that 
gives the most significant result and maximizes the drop of fRH 
between the two mass bins is Mum ~ 8 x 10 ^"^M©, for which 
fRH ^ 30% (40%) in the LM bin and fRH ^ 79% (79%) in the 
HM bin in the case (i) (in the case (ii)). For ^ Sx 10^^ Mq 
the observed fRH in the two mass bins differs from that obtained 
by the Monte Carlo analysis at ~ 3.2cr in the case (i) and ~ 2.5cr 
in the case (ii). This means that the chance probability of the 
observed drop of fRH is < 7.4 x 10“"^ (i) and < 5.7 x 10“^ (ii). 
Based on this analysis we conclude that there is statistical ev¬ 
idence for a drop of the fraction of RHs in galaxy clusters at 
smaller masses. A similar conclusion was obtained using X-ray 
selected clusters (Cassano et al. 2008), however this is the first 
time that such indication is derived using a mass-selected sam¬ 
ple. In Sect. 7 we discuss possible biases due to incompleteness 
in our current sample. 

6. Radio Halo-cluster merger connection 

In this Section we investigate the connection between the pres¬ 
ence/absence of RHs in clusters and the cluster dynamical status 
merger/relaxed. Following Cassano et al. (2010), in Fig. we 
report the cluster morphological parameters, derived in Sect. 
for the 50 clusters of the sample with available Chandra data, in 
three diagrams: c-w, C-P3 jP^, W-P3 jP^. Vertical and horizon¬ 
tal dashed lines are taken from Cassano et al. (2010), these lines 
represented the median value of each parameter, and were used 
to separate merging (w > 0.012, c < 0.2 and F 3 /F 0 > 1.2x 10“^) 
and relaxed (w < 0.012, c > 0.2 and F 3 /F 0 < 1-2 x 10“^) clus¬ 
ters. Here we use these lines as a reference to compare our mea¬ 
surements with previous published result^ Fig.j^shows that RH 
clusters (red dots) can be separated from clusters without RH 
(black dots) in the morphological diagrams: RHs are associated 

^ We note that Cassano et al. (2010) derived these lines on a smaller 
redshift range (0.2 - 0.35), however there is no clear indication about an 
evolution of P3/P0 and w with z (WeiBmann et al. 2013b and references 
therein) 
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case No. 

Allim 

(10‘^Mo) 

Nrh 

(HM) 

^cluster 

(HM) 

Nrh 

(LM) 

^duster 

(LM) 

/rh 

(HM) 

Irh 

(LM) 

fti 

(HM) 

Z 

(HM) 

(i) 

6 

24 

50 

0 

7 

48% 

0 % 

42% 

2.49 

7 

16 

28 

8 

29 

57% 

28% 

42% 

2.29 

8 

11 

14 

13 

43 

79% 

30% 

42% 

3.17 

9 

5 

6 

19 

51 

83% 

37% 

41% 

2.07 

10 

4 

4 

20 

53 

100 % 

38% 

40% 

2.21 

(ii) 

6 

27 

50 

1 

7 

54% 

14% 

49% 

2.10 

7 

17 

28 

11 

29 

60% 

38% 

49% 

1.76 

8 

11 

14 

17 

43 

79% 

40% 

49% 

2.50 

9 

5 

6 

23 

51 

83% 

45% 

49% 

1.72 

10 

4 

4 

24 

53 

100 % 

45% 

48% 

1.93 


Table 2. The four clusters with suspect diffuse radio emission are considered as non RH clusters in the upper panel (/) and as RH clusters in the 
lower panel (//). 




[10^^ MJ [1014 MJ 

Fig. 5. Observed fraction of RHs fRH (dots) compared to the value predicted by the Monte Carlo simulations (shadowed regions) in the HM bin 
(left panel) and in the LM bin (right panel) as a function of the limiting mass In both panels the green dots and the green shadowed regions 
represent the case in which the 4 clusters with suspect diffuse emission are considered as RH clusters. 


with dynamically disturbed clusters, while the greatest majority 
of clusters without Mpc-scale diffuse radio emission are relaxed 
objects. About 80% of the clusters in the HM bin of our sample 
are mergers, and this explains why RHs are fairly common in 
this bin (Sect. |^. The only RH cluster that always falls in the 
region of relaxed clusters is A1689, however this cluster is un¬ 
dergoing a merger event at a very small angle with the line of 
sight {e.g., Anders son & Madejski, 2004), therefore its morpho¬ 
logical parameters are likely biased due to projection effects. We 
note that also clusters with relics and without RHs (blue dots) 
are located in the regions of dynamically disturbed systems, in 
line with literature observations {e.g., de Gasperin et al. 2014). 

We also note that at least 10 merging clusters of our sample 
do not host RHs. The existence of massive and merging sys¬ 
tems without RHs is well known (Cassano et al. 2010, 2013, 
Russell et al. 2011). If RHs are due to turbulence acceleration of 
relativistic electrons during cluster mergers they should have a 
typical lifetime of ~ Gyr (see Brunetti et al. 2009), which is of 
the same order of the merger timescale. However the generation 
(and cascading from large to smaller scales) of turbulence and 
its dissipation take some time, corresponding to a “switch-on” 
and “switch-off” phases that span a substantial fraction of a Gyr. 
This produces a partial “decoupling” between X-rays and radio 


properties, as during these phases RHs would appear underlumi- 
nous/absent whereas the hosting cluster would appear disturbed 
in the X-rays {e.g., Donnert et al. 2013). 

An additional possibility is that some of the dynamically dis¬ 
turbed systems host RHs with very steep spectrum, that are not 
easily seen at our observing frequencies (Cassano et al. 2006, 
Brunetti et al. 2008). Indeed the great majority of merging clus¬ 
ters without RHs belong to the LM bin, which might support the 
idea that in these cases (or some of them) the energy provided by 
the merger is not sufficient to generate RHs emitting at the ob¬ 
serving frequencies. In fact, this second possibility is expected 
to contribute to the drop of the fraction of RHs in less massive 
systems (Cassano et al. 2010, 2012), as currently observed in our 
sample (Sect.j^. 

7. Sample completeness 

As explained in Sec. in the calculation of the occurrence of 
RHs we included only clusters with available radio information 
about the presence/absence of RHs: these are 21/21 clusters in 
the low-z sample (z < 0.2) and 36/54 in the high-z (z > 0.2) 
sample. Considering that the PSZ sample gives a completeness 
of ~ 90% for the low-z sample and 80% for the high-z sample. 
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Fig. 6. (a) c - w, (b) c - P3IP0, (c) w - P3/P0 diagrams. Vertical and horizontal dashed lines: c = 0.2, w = 0.012 and P3/P0 = 1.2 x 10"^. Red, 
black and blue dots represent clusters with RH, clusters without RH and clusters hosting relics (without RHs), respectively. Black open dots are 
clusters with suspect diffuse emission from the NVSS. 


we can estimate a completeness of our sample (which takes into 
account both the completeness in mass and in the radio informa¬ 
tion) of ~ 63^ 18 clusters in the high redshift range (z > 0.2) 
still lack radio information: 17 with M < 8 X 10^"^Mo and one 
with M > 8 X lO^^Mo. 

In this Section we evaluate how much the omission of these 
clusters can affect our results. 

We consider the total sample of 75 clusters and assume three 
extreme cases: 

a) all the missing clusters with M < 8 x 10 ^"^Mq are clusters 
without RH and the only one with M > 8 X IO'^Mq is a RH 
cluster; 

b) all the missing clusters with M < 8x10^4^0 host a RH, and 
the one with M > 8 x 10 ^"^Mq is a non RH cluster. 


^ This is estimated as (2i/o.9)+(54/o.8) ~ 63% 


c) the fraction of RHs in the 18 missing clusters is independent 

of the cluster mass. 

Cases (a) and (b) are simply adopted to obtain the maxi¬ 
mum (a) and the minimum (b) drop of the RH fraction with 
mass that can be expected starting from current data. We stress 
however that the case (b) is particularly unlikely since it implies 
that the occurrence of RHs is stronger in less massive systems, 
which is not justified by any observational results achieved so 
far (e.g., Cassano et al 2008). 

We examined the scenarios listed above in both the cases 
where the 4 low-z clusters with suspect diffuse emission are con¬ 
sidered as (/) non RH clusters and as (ii) RH clusters. 

In case (a) we add 17 non RH clusters to the LM bin and one 
RH cluster to the HM bin, thus the fraction of clusters with RHs 
in the LM bin becomes /rh = 13/60 = 22% (17/60 = 28%), 
while in the HM bin//?// = 12/15 = 80% (12/15 = 80%) in case 



















































































V. Cuciti et al.: On the occurrence of Radio Halos in galaxy clusters 


(i) (in case (ii)). Adopting the Monte Carlo approach, described 
in Sect.j^ we find that this corresponds to a 4.2 cr (3Jcr) result 
in case (i) (in case (ii)). 

In case (b) we add 17 RH clusters to the LM bin and one non 
RH cluster to the HM bin, so that fRuiLM) = 30/60 = 50% 
(34/60 = 57%) 2in&fRH{HM) = 11/15 = 73% (11/15 = 73%), 
with a 1.6 cr (1.2 cr) significance level in case (i) (in case (ii)). 

Finally, in case (c), we assumed that the fraction of RHs in 
the 18 missing clusters is the same we measured in the sample 
of 57 clusters analysed in the present paper: ~ 42% (/), ~ 49% 

(ii) , independently of the cluster mass. The fraction of cluster 
with RHs in the LM bin would be fRuiLM) = 20/60 = 33% 
(25/60 = 42%) and in the HM bin fRH(HM) = 11/15 = 73% 
(11/15 = 73%), corresponding to a 2.8 cr (2.2cr) result in case 
(i) (in case (ii)). 

Based on our analysis we conclude that the evidence for a 
drop of the fraction of clusters with RH at smaller masses is 
tempting and cannot be completely driven by possible biases de¬ 
riving from the incompleteness in mass of the (radio) sample. 
Namely, even in the very unlikely and extreme case (b), a hint 
of difference in the occurrence of RH still remains between the 
high-mass and low-mass systems in our sample. 

8. Summary & Conclusions 

The study of the statistical properties of RHs in galaxy clusters 
has became increasingly important in the last decade: it is a pow¬ 
erful tool to test the theoretical models for their origin and to un¬ 
veil the connection between RHs and cluster formation. In their 
simplest form, homogeneous re-acceleration models predict that 
RHs should be found in massive and merging objects, whereas 
the fraction of clusters with RHs, /rh, should drop towards 
smaller merging systems and RHs should be absent in relaxed 
clusters. In order to test these expectations large mass-selected 
samples of galaxy clusters are necessary. SZ-cluster surveys, i.e. 
the Planck-SZ survey (PSZ, Planck Collaboration 2014) have re¬ 
cently enabled the construction of cluster samples that are almost 
mass-selected, thanks to the tight correlation between the SZ sig¬ 
nal and the cluster mass (Motl et al. 2005; Nagai 2006). Recent 
studies, based on the EGRHS (Venturi et al. 2007,2008, Kale et 
al. 2013,2015) and the PSZ catalogue, have shown the presence 
of a bimodal split between clusters with and without RH also in 
the radio-SZ diagrams for clusters with Tsoo > 6 x 10“^ Mpc^ 
(Cassano et al. 2013). However, the mass completeness of the 
sample used by Cassano et al. (2013) is 50% and does not allow 
to probe the existence of a drop in /rh towards small clusters. 

Here we have presented a step toward an unbiased analy¬ 
sis of the occurrence of RHs, as a function of the cluster mass, 
in a mass-selected sample of galaxy clusters. We built a sam¬ 
ple of 75 clusters with M > 6 x 10in the redshift range 
0.08 <z< 0.33 selected from the Planck SZ catalogue. Among 
these clusters 57 have available radio information, for 21/21 
clusters in the redshift range 0.08 - 0.2 we used NVSS and liter¬ 
ature information, whereas 36/54 clusters at z = 0.2 - 0.33 have 
data from the EGRHS (plus literature information). Our study 
is based on these 57 clusters. The completeness in mass of this 
sample is ~ 63%, larger than that available in previous studies 
(e.g. Cassano et al. 2013). We also used the available Chandra 
X-ray data, for 50 out of 57 clusters, to derive information on the 
cluster dynamical status. 

The presence/absence of RHs has been determined by using lit¬ 
erature information for all the high-z clusters (z > 0.2, Venturi 
et al. 2007,2008; Kale et al. 2013,2015 and references therein), 
and for the majority (14) of low-z clusters (z < 0.2). We repro¬ 


cessed and analysed NVSS data of the remaining 7 low-z clus¬ 
ters that lack literature information and conclude for possible 
diffuse emission in 4 cases (Sect. |3.l| ). We split our sample into 
two mass bins, the low mass bin (LM, M < Mum) and the high 
mass bin (HM, M > Mum) and derived the fraction of clusters 
with RHs in the two mass bins for different values of M//^, find¬ 
ing that /rh is 60 - 80% in the HM bin and 20 - 30% in the 
LM one. By means of Monte Carlo simulations we obtained the 
distributions of RHs in the two mass bins (after 10^ trials), ex¬ 
pected in the case that RHs were distributed independently of the 
cluster mass. We found that for Mum ^ 8 x 10the observed 
/rh in the two mass bins differs from that obtained by the Monte 
Carlo analysis with a significance that ranges between 2.5 cr and 
3.2cr, which means that the probability to obtain by chance the 
observed drop of /rh is < 5.7 x 10"^ or even lower (see Sect[^. 
This highlights the statistical significance of our results and sug¬ 
gests that the increase of the occurrence of RHs with the cluster 
mass is likely to be real, rather than by chance. 

The possibility of a drop of the fraction of clusters hosting 
RH for less massive systems is particularly intriguing. Indeed 
this is naturally and uniquely expected in the framework of tur¬ 
bulent re-acceleration models (e.g., Cassano & Brunetti 2005) 
that provide a popular picture for the formation of giant RHs in 
galaxy clusters. A solid comparison between models and our ob¬ 
servations is still premature due to the incompleteness of the ob¬ 
served sample (Sect.j^. Still, with this caveat in mind, in Fig.|^ 
we compare our measurements of the occurrence of RHs in the 
two mass bins (black solid line) with the formation probabil¬ 
ity of RHs derived from the turbulent re-acceleration model in 
its simplest form (red line). Specifically, following Cassano & 
Brunetti (2005), we adopted the semi-analytic Press & Schechter 
theory (PS, Press & Schechter, 1974) to generate merger-trees 
and follow the hierarchical evolution of galaxy clusters through 
merger events. We assumed that a fraction, rjt of the PdV work 
done by the infalling subclusters during mergers is channelled 
into magnetosonic waves that accelerate relativistic electrons, 
which in turn emit synchrotron radiation. We calculated the the¬ 
oretical evolution of /rh with the cluster mass in the redshift 
range z = 0.08 - 0.33 for given values of the model parame¬ 
ters (see caption of Fig. |7] for details). Uncertainties on the pre¬ 
dicted formation probability are estimated by running MC ex¬ 
tractions from the pool of theoretical merger trees and account¬ 
ing for the statistical variations that are induced by the limited 
size of the two observed subsamples defined in Tab. (using 
Miim = 8 X 10^"^ Mq). Despite the crude approximations adopted 
in these models, there is an overall agreement between the ob¬ 
served and predicted behaviour of fRu with the cluster mass. The 
model slightly underestimates fRu in the high mass bin. This 
may be due to two main reasons: 1) the use of the PS formal¬ 
ism, which is well know to underestimate the merging rate, and 
hence the number density, of very massive systems; 2) the fact 
that the model predictions do not include RHs with very steep ra¬ 
dio spectra, i.e., those with steepening frequency < 600 MHz. 
As an example in Fig.|7](black dashed lines) we report the effect 
on the observed statistics induced by removing from the sample 
USSRHs (marked with in Tab.[T]) and candidate USSRHs 
(marked with “c” in Tab. [T]), for which we don’t know the de¬ 
tailed spectral shape. 

With the procedure described in Sec. we analysed the 
Chandra X-ray data of 26 clusters and we derived the morpho¬ 
logical parameters (the centroid shift, w, the power ratio, P 3 /P 0 
and the concentration parameter, c), which are powerful diag¬ 
nostic of the cluster dynamical status. We combined them with 
previously published results (Cassano et al. 2010, 2013) and we 
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Fig. 7. Expected fraction of clusters with RHs with steepening fre¬ 
quency (Cassano et ah, 2010) > 600 MHz in the redshift range 

0.08 < z < 0.33 (red line and shadowed region). Calculations have 
been performed for the following choice of model parameters: b - 1.5, 
{B) -1.9 jdG (where B - (B) x (M/(M))^) and rjt = 0.2 (see Cassano et 
al. 2012 and referencees therein). The observed fraction of clusters with 
RHs in the two mass bins is also shown (black points with horizontal 
errorbars). The balck points with dashed errorbars show /rh in the case 
we exclude USSRH (and candidate USSRH) from the observed sample. 


confirmed that RHs are hosted by merging clusters, while the 
majority of non-RH clusters are relaxed, thus highlighting the 
key role that merger events play in the origin of RHs. We note 
the presence of few merging clusters without RHs. This obser¬ 
vational fact adds constraints for the origin and evolution of RHs 
that have been briefly discussed in Sect. in the context of cur¬ 
rent models. 

The calculation of the occurrence of RHs has been per¬ 
formed only for clusters with radio information about the pres¬ 
ence of diffuse radio emission in form of RH. 18 clusters are not 
included in our analysis because they still lack radio informa¬ 
tion. In Sec. 1^ we tested the possible effects of the sample in¬ 
completeness on our results assuming three different situations 
for the derivation of the final /rh (see Sect. [7] for details). We 
found that even in the most unfavourable case, that however is 
very unlikely (case (b) in Sect. [7]), a drop of the fraction of RHs 
at smaller masses would still remain. 

This is the first step of this study; observations of the missing 
clusters with the GMRT and the VLA are already in progress and 
will allow the conclusive measure of the occurrence of RHs in a 
mass-selected sample of galaxy clusters. 
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